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We investigate observational constraints on the generalized Chaplygin gas (GCG) model as the 
unification of dark matter and dark energy from the latest observational data: the Union SNe la 
data, the observational Hubble data, the SDSS baryon acoustic peak and the five-year WMAP shift 
parameter. It is obtained that the best fit values of the GCG model parameters with their confidence 
Os| ! level are A s = 0.731^6 ( 1<T ) ±aog ( 2<T )> a = -°- 09 ±o.}2 ( lcr ) -oil ( 2cr )- Furthermore in this model, 

we can see that the evolution of equation of state (EOS) for dark energy is similar to quiessence, 

Ph, 

and its current best-fit value is wode = —0.96 with the lcr confidence level —0.91 > wode > —1.00. 
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1. Introduction 

6 : n 

The recently cosmic observations from the type la supernovae (SNe la) the cosmic microwave background 
(CMB) Q], the clusters of galaxies Q etc., all suggest that the expansion of present universe is speeding up rather 
than slowing down. And it indicates that baryon matter component is about 5% for total energy density, and about 
95% energy density in universe is invisible. Considering the four-dimensional standard cosmology, the accelerated 
expansion of the present universe is usually attributed to the fact that dark energy (DE) is an exotic component 
C*~) ■ with negative pressure. And it is shown that DE takes up about two-thirds of the total energy density from cosmic 
observations. Many kinds of DE models have already been constructed such as ACDM quintessence Q], phantom 
quintom Q, generalized Chaplygin gas (GCG) modified Chaplygin gas holographic dark energy (icj . 



agegraphic dark energy ll| , and so forth. Furthermore, model-independent method 1 and modified gravity theories 
(such as scalar-tensor cosmology [ljj], braneworld models [3|) to interpret accelerating universe have also been 
discussed. 



■ ■ - 0. 



It is well known that the GCG model have been widely studied for interpreting the accelerating universe 
The most interesting property for this scenario is that, two unknown dark sections in universe-dark energy and dark 
matter can be unified by using an exotic equation of state (EOS). In this paper, we use the latest observational data: 
the Union SNe la data j^], the observational Hubble data (OHD) [21]], the baryon acoustic oscillation (BAO) peak 



from Sloan Digital Sky Survey (SDSS) [2g and the five-year WMAP CMB shift parameter [23| to constrain the GCG 



'Electronic address: lvjianbo819@163.com 

1 Using mathematical fundament, one expands equation of state of DE Wd e or deceleration parameter q with respect to scale factor a or 
redshit z. For example, Wd e (z) = i«o=const |12|. Wd e (z) = + wi z ll3l . Wd e ( z ) = w + w i l n (l + z ) H3 . w<je(- 2 ) = w + f^f Il5|l . 
q(z) = go + H z [l2l| . q{z) = qo + j¥j Q i where wo, itfi, or 50, Qi are model parameters. 
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model. And we discuss whether the parameter degeneration 24] 25] for the GCG model can be broken by the latest 
observed data, since it is always expected that the model degeneration problem can be solved by the more accurate 
observational data. 

The paper is organized as follows. In section 2, the GCG model as the unification of dark matter and dark energy 
is introduced briefly. Based on the observational data, we constrain the GCG model parameter in section 3. The 
evolutions of EOS of DE and deceleration parameter for GCG model are presented in section 4. Section 5 is the 
conclusions. 



2. generalized Chaplygin gas model 

The GCG background fluid with its energy density pgcg and pressure pgcg are related by the EOS 

v = -j a , (1) 

where A and a are parameters in the model. When a = 1, it is reduced to the CG scenario. 

Considering the Friedmann-Robertson- Walker (FRW) cosmology, by using the energy conservation equation: 
d(pa 3 ) — —pd(a 3 ), the energy density of GCG can be derived as 



Pgcg — Pogcg 



[A s + (l-A s )(l + zf^}Th, (2) 



where a is the scale factor, A s — {+ a . For the GCG model, as a scenario of the unification of dark matter and 

Po 

dark energy, the GCG fluid is decomposed into two components: the dark energy component and the dark matter 
component, i.e., pgcg = Pde + Pdm, Pgcg = Pde- Then according to the general recognition about dark matter 

Pdm = P0dm(l + z) 3 , (3) 

the energy density of the DE in the GCG model is given by 

Pde = PGCG - Pdm = P0GCg[A s + (1 - A s )(l + z) 3 ( 1+Q )]TT^ - p adm (l + Z ) 3 . (4) 

Next, we assume the universe is filled with two components, one is the GCG component, and the other is baryon 
matter component, ie., p t — pgcg + Pb- In a flat universe, making use of the Friedmann equation, the Hubble 
parameter H is expressed as 

H 2 = = H 2 E 2 = H 2 {(1 - n ob )[A s + (1 - A s )(l + zf^]rh + f> ob (l + z f}. (5) 

Where Hq = lOO/i km 5 _1 Mpc _1 is the present Hubble constant, h = 0.72 ± 0.08 is given by Hubble Space Telescope 



key projects 2fj. figb is the present value of dimensionless baryon matter density, and a joint analysis of five-year 



WMPA, SNe la and BAO data gives fl 0b h 2 = 0.02265 ± 0.00059 
observations to constrain the GCG model parameter ( A s , a) . 



271 ] . In the following section, we will use the cosmic 
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3. Constraint on GCG model parameter 

It is necessary for the investigation of type la supernovae to explore dark energy and constrain the models. Since 
SNe la behave as excellent standard candles, they can be used to directly measure the expansion rate of the universe 
up to high redshift with comparing with the present rate. Theoretical dark-energy model parameters are determined 
by minimizing the quantity [28[ 

2 fa s (Pobs(zt) - Vth{0; Zj)) 2 

XSNe(°) = Z^ ~2 ' ( 6 ) 

i—\ obs\i 



where N = 307 for the Union SNe la data [201 ]. which includes the SNe samples from the Supernova Legacy Survey 



32] and several 



[29|, ESSENCE Surveys [30j, distant SNe discovered by the Hubble Space Telescope [3l|, nearby SNe 
other, small data sets. The la error a bs;i are from flux uncertainty, intrinsic dispersion of SNe la absolute magnitude 
and peculiar velocity dispersion, which are assumed to be Gaussian and uncorrelated. 9 denotes the model parameters. 
fi b s is the observed value of distance modulus and can be given by the SNe dataset. The theoretical distance modulus 
fith is defined as 

(J-th(zi) = mth(zi) - M. (7) 
Here mth(z) is the apparent magnitude of the SNe at peak brightness 

m th (z)=M + 5log w (D L (z)), (8) 
and absolute magnitude M can be given by relating to the magnitude zero point offset M, 

M = M + (9) 

TT-l 

with (j,q = 5logio(j^) + 25 = 42.38 — blogioh. Thus according to Eqs. J7]), ([5]) and ©, the theoretical distance 
modulus can be written as 

IHh{z) = 5log w (D L (z)) + /i , (10) 
where Dl{z) is the Hubble free luminosity distance 

D L {z) = H d L {z) = (1 + z) f (11) 

Since the nuisance parameter /i is independent of the data and the dataset, from above equations one can see that 
the distance modulus of different SNe (i.e. at different redshift z), n(zi) and fJ-(zj) are uncorrelated. So, the covariance 
matrix included in the XsNe 0^1- ©) is diagonal with entries Cj. 

Furthermore, by ex pan ding the XsNe °f expression ([6]) relative to /xo, the minimization with respect to fj,o can be 
made trivially 28( 3 34 j 



X 2 SNe(0) = A(6) - 2^B{6) + filC, (12) 

where 

A (0) = ^ obs ^ - t^h^iiE° f °' 6 ')] 2 (13) 

i=l (T ' 



z 


0.09 


0.17 


0.27 


0.40 


0.88 


1.30 


1.43 


1.53 


1.75 


H(z) (kms" 1 Mpc)" 1 


69 


83 


70 


87 


117 


168 


177 


140 


202 


la uncertainty 


±12 


±8.3 


±14 


±17.4 


±23.4 


±13.4 


±14.2 


±14 


±40.4 



Table 1. The observational H(z) data [21 ] [37 1 



N 



B{0) = £ 



Vobs(zi) - fi t h{zi;no = 0,0) 



(14) 



C 



N 1 



(15) 



Evidently, Eq. ^ has a minimum for jiQ = B/C at 



X 2 SNe (9)=A(6)-B(9r/C. 



(16) 



Since XsNe min = XsNe min an d XsNe * s independent of nuisance parameter /x [34| , here we utilize expression ([To]) 
to displace ([6]) for SNe constraint. 

Since the Hubble parameter H(z) depends on the differential age of the universe, 

1 dz 



H(z) = 



1 



dt ' 



(17) 



the value of H(z) can directly be measured through a determination of dz/dt. By using the differential ages of 
passively evolving galaxies from the GDDS (sfj and archival data 36 1. Ref. [^l] got nine values of H(z) in the range 



of < z < 1.8 (see Table 1). Here the observed Hubble data H(zi) and H(zj) are uncorrelated, for they are obtained 
by the observations of galaxies at different redshift. Using these nine observational Hubble data one can constrain 
DE models by minimizing 



37] 38] 



XHub( H 0, 



y-y [H t h(Ho,0, Zi) — H bs(zi)] 2 



(18) 



1=1 



obs-i 



where H t h is the predicted value of the Hubble parameter, H b s is the observed value, a bs;i is the la uncertainty of the 
measurement of standard deviation. Here Hq contained in the xjj u b(^o^) as a nuisance parameter is marginalized 
by integrating the likelihood L(6) — J dHoP(Ho) exp (—x 2 {Ho,6)/2). P(Hq) is the prior distribution function of the 
present Hubble constant, and a Gaussian prior H — 72 ± 8 km S~ l Mpc -1 [2 61 ] is adopted in this paper. 

Using a joint analysis of Union SNe la data and OHD (i.e., Xtotai = XsNe + X% u b)> Fig. 1 shows the constraint 
on GCG parameter space A s -a at the \a (68.3%) and 2a (95.4%) confidence levels. For this analysis the best fit 
parameters are A s — 0.80 and a — 0.42. It is obvious that two model parameters, A s and a, are degenerate. And it 
can be seen that model parameter a has the lager variable range. Then in order to get the stringent constraint and 
diminish systematic uncertainties, in what follows we combine the standard ruler data (the BAO peak from SDSS 
and the five-year WMAP CMB shift parameter R) with the Union SNe la data and the OHD to constrain the GCG 
model. 

Because the universe has a fraction of baryon, the acoustic oscillations in the relativistic plasma would be imprinted 
onto the late-time power spectrum of the nonrelativistic matter 39j. Then the observations of acoustic signatures 
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Fig. 1. The 68.3% and 95.4% confidence level contours for A a versus a from the Union SNe data plus the OHD. 




in the large-scale clustering of galaxies are very important for constraining cosmological models. From the BAO 
constraint, the best fit values of parameters in the DE models can be determined by constructing 

[A{6) 



XbaoV) 



Where 



A{6) = ^/fk^E{z BAO r 1/3 [ 



1 



dz 



12/3 



zbao Jo E(z ;6) 

--n ob + {i-n ab )(i-A s )Th 



(19) 



(20) 



24] 25] 41]. The observed 



ilom is the effective matter density parameter given by ilo m = 

value A bs with its ler error oa is A Q b s = 0.469(n s /0.98) _a35 ± 0.017 measured from the SDSS at zbao — 0.35, here 



27] 



n s is the scalar spectral index [42j and its value is taken to be 0.96 as shown in Ref. 

The structure of the anisotropics of the cosmic microwave background radiation depends on two eras in cosmology, 
i.e., the last scattering era and today. They can also be applied to limit DE models by minimizing 43] 

(R(9) - R obs ) 2 



Xcmb(8) = 



Where the shift parameter 44 1 



R(6) = V^or, 



dz 



(21) 



(22) 



z rec — 1089 is the redshift of recombination. The observed value R bs = 1-710, and its corresponding lcr error is 
(7r = 0.019 according to the five-year WMAP result 0. 

Above four observational data are uncorrelated for each other, since they are given by different experiments and 



methods. Then the total likelihood 



Xtotal 



can be constructed as 



v 2 

Atotal 



XSNe 



XHub 



Xbao 



XCMB- 



/ 2 . 

v /;/ / /; 



(23) 
= 322.87, 



Using Eq. (|23]l we get the best fit values of GCG model parameters (A s ,a) are (0.73,-0.09) with \ 
and the reduced x 2 value is 2 Xmm/d°f = l-03. The lcr and 2a confidence level contours of GCG model parameters 



2 The value of dof (degrees of freedom) for the model equals the number of observational data points minus the number of parameters. 
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(a) (b) 
Fig. 2. The 68.3% and 95.4% confidence level contours for A s versus a from the Union SNe data plus the BAO data (a), and 

a combined analysis of the Union SNe, OHD, BAO and CMB data (b). 



are plotted in Fig. 2 (b). From this figure, we obtain the values of model parameters with the confidence levels, 



A s = 0.73l£^ (lcr) Iq 09 (2a) and a = -0mZ 012 (la) Z^ f 9 (2a). It can be seen that parameters A s and a are also 
degenerate, and at the la confidence level these results are consistent with the standard dark energy plus dark matter 
scenario (i.e., the case of a = 0). Furthermore, one can see that this constraint on parameter a is more stringent 
than the results in Refs. 



24] 25], where the constraint results for the GCG model parameters are A s = 0.70to \r and 



a = —0.09^^33 at 2a confidence level by using the X-ray gas mass fractions of galaxy clusters and the dimensionlcss 
coordinate distance of SNe la and FRIIb radio galaxies 24], and A s = 0.75*q'o8! a ~ ^-^^aVe a t %a confidence level 



by means of the observational Hubble data, the 115 SNLS SNe la data and the SDSS baryonic acoustic oscillations 
peak [25(. 

At last, we also consider the constraint on the GCG model parameter from a combination of Union SNe la and 
BAO data, the best fit happens at A s = 0.75 and a — 0, which can be reduced to the standard dark energy plus 
dark matter scenario. But at their confidence levels, the two parameters are also highly degenerate. In Fig. 2 (a), we 
display the constraint result for this analysis. 



4. Constraint on EOS of dark energy and deceleration parameter 

According the Eq. ((5]), deceleration parameter q in the GCG model can be obtained by 

si dH 

q = (l + Z )---l. (24) 

The equation of state of dark energy is derived as 

Pde -(l-n 06 K[A, + (l-^)(l + ;0 3 < 1+a >]-rf= 
Wde = = 1 ) (25) 

Pde (1 - Q ob )[A s + (1 - A s )(l + z)3(l+«)]— - fi Qdm (l + Z )3 



where f^odm is the present value of dimensionless dark matter density. Based on Eqs. (|24|) and (|25|) . the confidence 
levels of the best fit Wde(z) and q(z) calculated by using the covariance matrix are plotted in Fig. 3. From Fig. 3 
(a), it is easy to see that the best fit value wode = Wd e (z = 0) = —0.96 > —1, and the la confidence level of wode 
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Fig. 3. The best fits of Wde{z) and q(z) with ler confidence level for GCG model. 



is —0.91 > wode > —1.00. In addition, it can be found that the best fit evolution of u>d e (z) for GCG is similar to 
the quiessencc model (wd e (z) = const ^ —1). From Fig. 3 (b), we can see that the best fit values of transition 
redshift and current deceleration parameter with confidence levels are zt = 0.74+q'os (ler), qo = — 0.55^0 06 (Ic)- One 
knows that zt describes the expansion of universe from deceleration to acceleration, and qo indicates the expansion 
rhythm of current universe. Comparing our results with Ref. [45J], where zt = 0.491q'q7 (ler) and qo = —Q. 73^ 20 
(ler) are obtained from Union SNe la data by using a linear two-parameter expansion for the decelerating parameter, 
q(z) = qo + q\z, it is clear for our constraint, that the universe tends to an earlier time to acceleration and a milder 
expansion rhythm at present. 



5. Conclusion 

The constraints on the GCG model as the unification of dark matter and dark energy are studied in this paper by 
using the latest observational data: the Union SNe la data, the observational Hubble data, the SDSS baryon acoustic 
peak and the five-year WMAP shift parameter. We find that the model parameters A s and a are degenerate, and 
their values are constrained to A s = 0.73±°;gg (ler) torn ( 2a ) and a = -°- 09 ±ai2 ( lcr ) -oil ( 2cr )- This constraint 

24\ 3- 

Furthermore, it is shown that the evolution 



on parameter a is more stringent than the results in Refs. 
of EOS of dark energy for the GCG model is similar to quiessence, and the best fit value of current EOS of DE 
wode — —0.96 > —1. And it indicates that the values of transition redshift and current deceleration parameter are 
z T = 0.74t° ;E>5 (la), qo = -0.55±° .° 6 5 (la). 
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